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Abstract 

We find general solutions for the dimensionally regularised scalar one loop three-point and 
four-point functions with one heavy quark propagator. The scalar one-point function vanishes, 
. while the expression for the two-point function has been found before. For the latter we give a 

' detailed derivation. We also discuss some special cases and compile useful formulae. 

(N 

1 Introduction 

Since its early applications [Q the heavy quark symmetry has been one of the key ingredients 
in the theoretical investigations of hadrons containing a heavy quark. It has been successfully 
. applied to the heavy hadrons spectroscopy, to the inclusive as well as to a number of exclu- 

sive decays (for reviews of the heavy quark effective theory and related issues see or |^). 
^ ' To describe interactions with not too energetic light mesons, the heavy quark symmetry has 

Xjf^ \ been combined with chiral symmetries leading to the heavy hadron chiral perturbation theory 

m ■ (HChPT) 0. 

\ The one loop calculations within the heavy quark effective theory are considerably simplified 

■ if the light-quark masses are neglected. Very common in the HChPT is a similar approximation, 

with the finite contributions omitted, while only the leading logs are retained §,§. This 
' approximation is, however, becoming less adequate with the increasing precision of data coming 

from B factories. To go beyond the leading log approximation and/or take into account the 
^ • counterterms appearing at the next order in the chiral expansion, the general solutions for the 

one loop scalar functions need to be considered. In the context of the HChPT a general solution 
for the one loop scalar two-point function with one heavy quark propagator has been found in 
pH \ 07 Hi- We extend this calculation and find solutions for the scalar three-point and four-point 

functions with one heavy quark propagator. 

The vector and tensor one loop functions can then be expressed in terms of the scalar one 
loop functions using the algebraic reduction 13] . Also, the one loop scalar functions with two or 
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X 



\ more heavy quark propagators can be expressed in terms of the one loop scalar functions with 

just one heavy quark propagator. For the case of the equal heavy-quark velocities this can be 



accomplished using the relation 

1 1 



w-(7 — A u-g — A' A — A' \u-(7 — A v-q 



(1) 



For unequal heavy quark velocities techniques developed in [ pl)[ can be used. 

The scalar one-loop functions with heavy quark propagators can be derived also directly 
from the scalar functions of the full theory by using the threshold expansion (see also 
appendix B of ) . This technique has recently been used for the calculation of the scalar and 
tensor three-point functions with one and two heavy quark propagato rs p^,p4| . In the present 
paper, however, we will not follow the approach of Bouzas et al. [p0[,|13{,]14| but rather do the 
calculation from scratch. 

The paper is organized as follows: in section 2. we make some general remarks and list 
useful relations used further on in the calculation. In section 3. we calculate one and two point 
functions. We continue with the calculation of the three-point function in section 4. and with 
the four-point function in section 5. 
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2 General remarks and useful relations 



Before we start with the calculation, let us first list some useful relations and the conventions 
that are going to be used further on. The greater part of this section is a review of the relations 
and the conventions used in with certain modifications. The major difference between the 
conventional one loop scalar functions and the one loop scalar functions with one heavy quark 
propagator is the appearance of the propagator linear in the integration variable q. Therefore a 
modified version of the standard Feynman parameterization is used 

1 _.n n.. f^... ^(i-E."on.ew 



poo n j'' 

=Nl / 2dX Yl du. 



/o°° 2rfA /o dxi /(f ' dx2 • • • /(f dxN-i 



--N\ 



- xi) + Aj{xj^^ - X,) + Anxn-1 + 2SA]^+i 



where Q{u) is the Heaviside function 0(u) = 1 for it > and zero otherwise. In the calculation 
Ai are going to be "full" (inverse) propagators {{q + — mf + i5) and B the heavy quark 
propagator {vq— /S.+i6) . Note also, that the leading power of q^ in the denominator has increased 
from the left-hand side's (g^)^"*"^/^ to the right-hand side's {q^)'^'^^- The integration over q has 
been made more convergent, but then another integration over infinite range (integration over 
A) has been introduced through the parameterization. 
A very useful identity used in the calculation is 



a 



[{q + PiY - ml + i6][{q + P2Y - ml + 16] [[q+piY - ml + i6][{q + lY - NP + i6] 



(3) 



[[q + P2Y -ml + iS] [[q + If - + i5] ' 



where a is an arbitrary parameter and 

l^pi+a{p2~pi), (4) 
Af2 = (1 - a)ml + ami " "(1 " ")(P2 - Pif ■ (5) 

The parameter a can then be chosen at will. It is useful to keep it real, though. Then there 
are no ambiguities connected with the shift of the integration variable that is performed, 
as usual, before the Wick rotation. For instance a can be chosen such that AP — 0. If 
{p2 — Pif < {mi — 7712)^ or {p2 — piY ^ (™i + ^2)^, then a is real. If one of the masses is 
made to be zero, the integration is simplified considerably (as will be seen in the calculation of 
the four-point function (|39|)). The other option used below is to set a such that P — 0. This 
can be done for real a if (but not only if) one of pi, p2 or pi ±p2 is timelike. This shows, that 
in general product of propagators at least one internal or one external mass can always be set 
to zero, even with a restricted to be real. 

In doing the integrals the following procedure proves to be very useful. Consider 

2d\ / dx—^ — ^ -. (6) 

The integration over x can be simplified by the change of the integration variables A = A' + 
where /3 is chosen such, that the coefficient in front of vanishes, i.e. /3 has to solve the 
equation 6/3^ + c(3 + a = 0. Then the integrand is linear in x, so the integration over x is trivial. 
The integration bounds are 

.1 j-i 
2dX I dx-- - ^ dx I 2dX' - - - 



Jo Jo J^0x , , 

2dX dx---+ 2dX dx---. 



Jo J-/3 J-\/l3 
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As the results of the mtegration, the functions such as logarithms, dilogarithms (Spence 
functions) and hypergeometric functions will appear. Since the arguments of the functions will 
in general lie in the complex plane it is necessary to discuss the conventions used. The logarithms 
used in this paper have a cut along the negative real axis. For x exactly negative real we use 
the prescription ln(x) ln(a; + ie), where e > is a positive infinitesimal parameter. In other 
words, ln(a;) = In \x\ + iir for x negative real^. In particular ln(— 1) is defined to be ln(— 1) = iir. 
Of course this choice is completely arbitrary and at the end of the calculation one has to check 
that results are independent of this choice. Using this definition for the logs of the negative real 
arguments the logarithm of an inverse is 

In = -ln(.x) + 27riRe(-^(x), (8) 

with 

^ I 1 ; a; on negative real axis, 

1 ; otherwise. 

Note that the change from the usual rule for the logarithm of an inverse is just on the negative 
real axis. For the arguments away from the negative real axis the function Re*-^-* {x) is exactly 
zero and everything is as usual. The logarithm of a product is 

ln(a6) = ln(a) + ln(6) +7?(a,5), (10) 

where 77 function is ^ 

27ri { 8 (- Im a) e (- Im 6) e (Im a6) 



?7(a,6) = < 



, , a and b not negative real, 
- e(Ima)e(Im6)e(-Ima6)} 

— 27ri{0(Ima) + 8(Im6)} ; either a or 6 negative real, 

— 27ri ; a and b negative real. 



(11) 



The normal rule for the logarithm of a product applies for these important cases 

ln(a6) = ln(a) + ln(6); Ima and Im& have the opposite sign, 
ln(a/6) = ln(a) — ln(6); Ima and Im& have the same sign, 

with a, b not negative real. 



(12) 



3 One- and two-point functions 

In this section we will concentrate on the calculation of the dimensionally regularized one-point 
and two-point functions in the heavy quark effective theory 

"16^^°^^^ " (2^ J \vq-A + t6y 

"le^^"^""' ^ " (2^ J '^{vq- A + i5){q^ -m? +i5y 

^Note that this prescription does not change tlie calculation of the logarithms away from the negative real axis. 
In particular it does not change the value of a logarithm with an argument that already has an infinitesimal but 
nonzero imaginary part. For more discussion on this point see text after Eq. (^). 

^Note, that in comparison with |l5[ , the 77 function has been extended also to the negative real arguments (cf. 
discussion after Eq. and Eq. (|68|)). For arguments away from the negative real exis (also if by an infinitesimal 
amount) it is the same as in |15|]. 



(13) 
(14) 
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where 5 is a positive infinitesimal parameter, n = 4 — e, while m and A are real. The general 
solution for the two-point function has been found by Stewart in Ref. |^ (see also |^ and 
references therein). In this section we will derive Stewart's result. 
We start with the integral 



Ir = 



(27r)" 

Using the Feynman parameterization 



{vq- A + iS){q'^ - m2 + iSy 



1 



dX- 



we get 



/, = fi' 



r(r)r(s) Jo (a + 25A)■'■^ 
(-1)-'- 2r(r+ 1 - f) 



r(r) 



(47r)^ 



dA(A^ + 2XA + m-' - iS) 2 



(15) 



(16) 



(17) 



The integral can be expressed in terms of the hypergeometric function. Introducing the new 
variable A' = A + A and then splitting the integration interval for negative A we get 



dX' 



(A'" - A2 



i6) 



N 



dX 



(A2 - A2 



iS) 



N 



e(-A) 



dX 



IA| (A2 - A2 + m2 - i,5) 
(18) 



N ■ 



where we write N 
{rn? - 



iS) 



n/2 for short. Another change of variables u 
du \ u T 



X leads to 



-N 



A2 
2 



-vn? — A2 
du r 



1 







2V^ 



-N 



e( 



-A)}. 



(19) 



These integrals can be expressed in terms of the hypergeometric functions 2F1 (a, /3; 7; z) through 
the identities ( 64a ), (341:) listed in the appendix Using the transformation formula ( |63|) 
together with 2-^1(0, (3; 7; z) — 2i^i(a, 0; 7; z) = 1 we arrive at 



Ir = /i' 



(-1)-'' 2r(A^) 
r(0 (47r)T 



- A{m^ - A? - i5) 
r (N 



-N 



FA N. 



(2N - i)r(iv) 



(m 



1 3 _ 

2' 2' m 

-i(5) 



-A2 



N 



(20) 



where N = r + 1 — n/2. 

Let us first discuss the case when r is equal to zero or negative integer, i.e. when the 
integrand, apart from the heavy quark propagator, is a polynomial. The integrals for the 
physical case n —^ 4 are divergent, but we can make sense of it through analytic continuation. 
At fixed r the integral Ir is taken to be an analytic function of the complex dimension n. For 
r equal to zero or negative integer and n/2 ^ Z all functions appearing in ( pO| ) are finite, apart 
from r(r) that is infinitely large. Thus Ir vanishes for r zero or negative integer everywhere in 
the n complex plane apart from the points on the real axis with integer n/2 ^ r + 1. Analytic 
continuation of Ir (j2^) is then equal to zero in the whole n plane. Integrals over polynomials 
(and one heavy quark propagator) are in the dimensional regularisation thus equal to zero. In 
particular, the one-point scalar function Aq = 0. 

For the two-point function we have r = 1 and therefore N = e/2. So the two point function 

is 



7r(m 



{-A){m^ - A^ -i6y 
'■-A^-^5)iy 



2F1 



1 1.3. -A2 

2'2'2'm2-A2-iJ 



(21) 
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Since r(e/2) 2/e — 7 + 0{e) we have to expand hypergeometric function around e/2 = in 
order to get the finite terms correctly 



e 1 3 



2' 2' 2' 



d 



2dN 



1 3 



2' 2' 



JV=0 



The partial derivative can be found using the series expansion m^) and is 



d 

dN 



'2'2' 



= -ln(l - z) 



N=0 



(22) 



(23) 



This leads us to the final result for the two-point function 



(27r 



,4-e 



1 



(v-q - A 
2A (2 



iS) 



I 7 + In 47r — In 



2-2F 



(24) 



where F(x) is a function as defined in 
always taken to be positive real) 



valid for both positive and negative A (while m is 



— \/ — \ ln(.x + \J x"^ — \ + i5) 

X 



TT 

tan 

2 



\x\ > 1, 
\x\ < 1, 



(25) 



with S an infinitesimal positive parameter. Note that for a; < — 1, F(l/x) has an imaginary part 
that corresponds to the particle creation. Also, the two point function has to be a continuous 
function of A, as can be seen from ( p^ or (p^). It is easy to check that for |x| = 1 the two-point 
function is continuous as then F{±1) ~ 0. The two-point function is also continuous for A ^ 0. 
Even though F{l/x) diverges as a; 0, the two point function ( p^ is finite and equal to m/SiT. 

Finally, for r ^ 2 both r(r) and T{N) in (|^) are finite in the hmit n ^ 4, so that Eq. (|^ 
can be used directly, with n set to n = 4. 



4 Three-point scalar function 

The one loop scalar three-point function with one heavy quark propagator is given by 

1 =, , , . N /" m 1 

;Co[v,k,A,mi,m2) = -— — / d q- 



167r2 ' ' ' " " i2n)" J ^{vq-A + iS){q' -ml +i6){{q + k)^^ml+ 15)' 

(26) 

This integral is finite in 4 dimensions, so that e can be set to zero. Using the Feynman param- 
eterization (||) we get 

~ Jo '^'^^ Jo [A2 + fc2x2 + 2vkxX - (fc2 +ml- ml^x + 2AA + - iS\ ' ^'^^^ 

The integration over x can be made trivial through the change of variables A = A' + ax. We 
choose a to be the solution of 

[k + avf=Q, (28) 



as then the term quadratic in x is zero. The solution is ai,2 — {—v-k± \J (u-fc)^ — fc^) and is real 
for any real four-vector k^ (this can be easily seen by going in the frame, where — (1, 0, 0, 0) 
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with the square root then equal to V fc^). Changing the integration order as in (0), and then 
integrating over a;, we get 



2dX ^^^fX^ + CX + D + {A\ + B] 



+ 



(AX + B) 
^ 2adX 



X^ + CX + D 

' / . X ^AHa^X^ -CaX + D + {-AaX + B)] 
{—AaX + B) I 

- \n[a^X^ - CaX + D + {-AaX + B)X] | 



A = 2{v-k + a), B ^2Aa + ml-ml- k'^, 
C = 2A, D = ml- iS. 



(29) 



where 



and a one of the sohitions a 1,2 of the quadratic equation (^8|). In ( |29| ) we have used the fact that 
D is the only complex parameter and split the logarithm in the second integrand. The integrands 
of both the first and the second integral have vanishing residua. In the second integral we then 
add and subtract the values of the logarithms for A = B/Aa and write 



2dX , fX^ + (A + C)X + (B + D) 
m ' 



X + B/A 

^ 2dX 
,X~B/{Aa) 



X^+CX + D 



In [a^X^ + C)aX + [B + D) 



-\ii[a'Xl~{A + C)aXo + {B + D)\) (30) 



2dA 



D\ 



- In [a{a - A)Xl + [B - Ca)Xo + £>] ) 



with Ao — B/{Aa). Note, that all three integrals in ( pO|) have integrands with vanishing residua. 
These integrals can be reduced into the sums of the dilogarithms, where care has to be taken 
regarding the imaginary parts of the arguments of the logarithms. The solutions of the integrals 
can be found in the appendix |^. The solution of the first integral can be found in (|5|), (|87|), 
with the definitions (^6|), ( |7^ ) (where ai = a2 — 0, note also the minus sign), while the solutions 
to the last two integrals can be found in (|69| ), ([7l|). Using the functions 5*3 and I2 defined in 
the appendix ^ the three-point function (^6[) finally reads 



Co 



2 r 
A 



h (0, Q,A + C,B + D,C, D, -B/A) 
+Sz (a^ - (A + C)a, B + D,B/Aa) 
-83 {a{a -A),B- Ca, D, B/Aa) 



(31) 



Note that the value of the three-point function in (pi] ) docs not depend on which of the solu- 
tions Q!i,2 of the equation ( p8| ) is used. This can be used as a useful check in the numerical 
implement at ion . 

The solution is simplified considerably if = 0. Then the x integration in (^7|) is trivial. 
Proceeding similarly as above we arrive at 



Co{v,k,A,mi,m2) 



v-k 



Li2 



Ao 



Ao — Ki 



■ln2(Ao 



(32) 
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where Aq — {m\ — m\)/{2v ■ k), while the solutions of 

+ 2{vk + A)A + 7712 -iS = {\- ki)(A - K3), 
A^ + 2AA + m^-i5 = (A-K2)(A-K4), ^^^^ 

and p{K,) = i-iy+\ 

The solution is even further simplified if besides = also nii = = ni- Then 



Co{v, fc, A, TO, to) = — 3- V p(Ki) 

1.2— n 1). h 



fe2=o v-k 2 



1 



(34) 



with Ki and p{ni) given in (|33| ). The three point function in this limit has been calculated before 
and is given explicitly in ]17|] (see Eq. (AlO) of [0). The two expressions agree completely. A 
number of numerical checks between numerically integrated expression ( p9| ) and final expression 
( |3l| ) have been performed as well. 

5 Four-point function 

The scalar four-point function with one heavy quark propagator is defined as 
1 



r£'o(w,Pl,P2, A,TOl,TO2,TO3) 



167r2 



(35) 



(2^)" J {vq ~ A){q^ - m?)((g +pi)2 - tuMq + P2? - ml) ' 



where the i5 prescription has been omitted in the notation. Again, the integral is convergent 
and e can be set to zero. Since the Feynman parameterization (^) is not symmetric in Ai and B, 
the elegant transformation used in the calculation of the conventional four-point function |l5|| 
and further improved in [p^ unfortunately cannot be applied. Instead, one repeatedly uses the 
propagator identity (^ to solve the integral (|35|). Since the parameter a in propagator identity 
(^) has to be real, the calculation differs depending on the values of the external momenta pi 
and p2 ■ 

First we take up the case, when one of the following inequalities is true Pi 2 ("^i + "^2,3)^ 
or 2 ^ ('^1 ~ ^2,3)^- If necessary, we renumber the momenta and reshuffle the propagators 
in (^Sj) in such a way that either p\ ^ (toi -t- TO2)^ or p\ ^ (toi — 777,2)^ in order to simplify the 
discussion. Then we use the propagator identity (^) on the second and the third propagators of 

(H) 



\-a 



((j2 -ml + i6){{q+pi)'^ - ™i + i5) [{q + Pi)^ -^1 + i6][{q + l)"^ - AP + i5] 

a 



(36) 



[g2 - ml + i5][{q + l)^ - AP + iS] ' 



where a is an arbitrary parameter and 

l^api, (37) 
AP = {I - a)ml + ami - a{l - a)pl. (38) 

We choose a such that — 0. This is satisfied by real a if either pi ^ (toi + TO2)^ or 
p1 ^ (toi — TO2)^ as has been assumed above. The scalar four point function is then 



(2^)4 J ^ (^;.g - A) [{q + p,)^ - m^] [{q + p,)^ - ml] [{q + l)^ + i6] 



« / ,4 a 
a q- 



(39) 



(2^)4 J ^ {vq - A)[q^ - TOf][(<7 + P2)2 - ml][{q + l)^ + i5] ' 
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with a the sohition of 

y{c? + {ra\ - mj ~ p\)a + m\^Q. (40) 
To calculate the two integrals in (^) it suffices to consider 

1 (41) 



(2^)4 J (v.g - A)[(9 + fci)2 - M2][((7 + k^Y - Mj][{q + fca)^ + iS] ' 

where the i5 prescription has not been written out explicitly in the first three propagators. 
Using the Feynman parameterization (j^) and integrating over q we arrive at 

/•oo /•! f^x 

Do= 2dX dx dy[p23X^ + pi2V^ + (pig - P23 - Pi2)xy + + (P2 - P3)a;A + P3A 

Jo "'0 "'0 

+ {Pi - P2)y\ + (-P23 + M|)a; + (p23 - Pi3 + Ml - Mi)y - iS]~^, 

(42) 



with 

Pij = {h - kj)'^, ^^^^ 



P, = 2{vh + A). 

To simplify the integration we introduce new variables y = xy' and A = xX' . The integration 
limits are then dx Jq dy 2dX — > dx xdy' x2dX' . Since x is positive and S an 
infinitesimal parameter of which only the sign matters, the extra factor of x'^ in the numera- 
tor can be canceled against the similar factor in the denominator. After the cancellation the 
denominator is linear in x. The integration over x is now trivial and yields 

/"OO /"l 

Do= 2dX / dy[{p23 - P13 + Ml ~ Ml)y + P3X - P23 + - *<5]"' x 
Jo Jo 

[Pi2z/ + A^ + (Pi - P2)yA + P2A + {-pi2 + Ml - Mi)y + M| - iS]-\ 

(44) 

To cancel the term in the integral above a new variable X = X' + f3y is introduced, with f3 
chosen to solve 

/32 + (Pi-P2)/3+pi2 = 0. (45) 
The solutions are real for (^2 — ki)^ real. We get 



1 /.O 

2 

13 



Do = { I 2dX ldy+ I 2dX / dy)- ^ (46) 

-A//3 [O'ly + oi A + cij [a2y + 02A + C2 + d2yX + A^J 



with 



ai = PP3 + P23 - Pi3 + Ml - Mi , a2 = [3P2 + Ml - M| - pu , 

b2=P2, 

C2 — M^ — iS, 

d2=2/3 + Pi-P2. 



bl=P3, h2=P2, 

ci = M| - P23 - iS, C2 = M2 - iS, 



(47) 
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After y integration we arrive at 



(Ai - A2) {bid2 - ai) 



+ 



2d\ 
A - Ai 
^ 2d\ 
A- A2 
2dA 
0A + A1//3 

2dA 

aTv^ 



In 



In 



biA + Ci 
61A + Ol + Cl 

feiA + ci 
61A + ai + ci 
(ai - /36i)A 



In 
In 



A^ + 62A + C2 



A2 



(62 + d2)A + 02 

A^ + &2A + C2 



C2 



Cl 



— /36iA + ai + Cl 

(«! — /36i)A + Cl 
—(3biX + ai + Cl 



A2 + (fe2+d2)A + a2+C2 

f3{P-d2)\^ + {a2-b2l3)X + C2 



(48) 



- In 



In 



/32A2-/3(62 + d2)A + a2 + c2 

/3(/3-d2)A2 + (a2-M)A + C2 
/32A2-/3(62 + rf2)A + a2 + C2 



with Ai_2 the solutions of 

(61^2 - ai)A^ + (02^1 - 0162 + cid2)X + a2Ci - aiC2 = (&id2 - ai)(A - Ai)(A - A2). (49) 

Note that the integrands above have vanishing residua, i.e. arguments of the two logarithms in 
the integrands are the same for A = Ai.2. Note as well, that the infinitesimal imaginary parts 
of Cl and C2, the —iS in (^7|), have to be equal. They originate from the same infinitesimal 
parameter in (^ ) that after the integration over x appears twice in (44). The size of the 
infinitesimal parts of A1.2 compared to —iS are thus unambiguously defined and have to be kept 
track of until the end of the calculation. 



The integrals 
in the appendix 



48| ) can be expressed in terms of the dilogarithms. This has been done 
. The solution for the first two integrals can be found in (^), with the 
definitions in ([76|), (|7^), while the solution for the last two integrals can be found in (fz^), 
d?!). Together with Eqs. (||), (|I|) and the cascade of abbreviations (|^), (||), (||) this gives 
the complete solution of the four point function with at least one external momentum pi, p2 
satisfying Pi 2 ^ ("^1 + "^2,3)^ or p2 ^ ^ (mi — 7712,3)^. Collecting the terms and rearranging 
the last two propagators in (^5|) if necessary, the four point function for pi ^ (mi + m2)^ or 
Pi ^ (mi — m2)^ finally reads 



Do{v,pi,p2,A,mi,m2,m3) =(1 - a)Do{v,pi,p2j, A,m2,'m3) 

+ a Do{v,0,p2,l,A,mi,m3), 

with I — api and a the solution of PiO^ + (m| — m^ — p^a + m^ — 0, while 

Do{v,h,k2,k3,A,Mi,M2) = \ w, , ^ r 

(Ai - A2) (oirf2 - ai) L 

I2 (-C1/61, -(ai + Cl)/5i, 62, C2, &2 + C?2, 0.2 + C2, Xl) 

-^2(-ci/6i, -(fli + ci)/bi,b2,C2,b2 + d2,a2 + C2, h) 

-h{0^ ~ I3d2,l3^,a2 - /352, -/3(62 + rf2),ai - /35i, -/35i, C2, 02 + C2,ci,ai + ci, 
+h{0^ - f3d2,P'^,a2 - (3b2, -f3{b2 + d2),ai - Pbi,-(3bi,C2,a2 + C2,ci,ai + ci, 



(50) 



(51) 



-A2//3)' 



with fli . . . 1^2 defined in (^7|), /3 defined in (^), p^j, defined in 



and Ai^2 solutions of (p^). 
and (E5h that in general 



Note that parameters a and /? are solutions of quadratic equations 
have two real solutions each. The value of the four point function Dq{v, ki, k2; k^, A, Mi, M2) 
does not depend on which of the solutions are chosen in the evaluation of (|50|) , (^l|) . This fact 
can be used as a useful test in the numerical implementation of the expressions given above. 

Now we take up the special case of pi on the light-cone, i.e. Pi = 0. From ( |40| ) it follows 
that a — ml/ {nil ~ ^2)- If "^1 "^2 then a is finite and the calculation proceeds as before. 
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(p5|)-(pq). For equal masses mi and m2, we evaluate the integrals by first taking mf ^ rrv^ and 
then performing the limit m\ — > m\ and thus \a\ — > oo. The external momenta in the last 
propagators of ( ^9| ) arc both equal to I — api . Thus the last external momentum in Dq of ( ^l| ) 
is going to he = I — api for both of the integrals in (^9|). In the limit a ^ oo the following 
leading order values ( ^ ) are obtained: ai — > 2api-{j3v — p2), hi 2av-pi, a 2api-p2 , where 
we have used also the fact that Pi = 0. The other coefficients 02, ^2, C2, d2 do not depend on a. 
The first term in the denominator of the integrand in ( |46| ) is then proportional to a, while the 
second term in the denominator does not depend on a at all. The a in the denominator cancels 
against the a in the numerator of (p^). For the case of equal masses = rn^ and p\ = 0, the 
solution is then the same as for the case 7^ "^2 except that (i) one has to replace (1 — a) and 
a in ( |39| ) with —1 and 1 respectively, and that (ii) oi, 61, ci in (|47|), ( p8| ) have to be replaced 
by their limiting values (divided by a) 

ai a[ = 2pi-{/3v - P2), 

h ^ b[ = 2vp,, (52) 
Cl ^ 4 = 2pi ■p2 - iS, 

where /? is the solution to (^^ (with Pi defined in (^)), and is different for the two integrals 
in (^). In the limiting value of ci coefficient given in ( |5^ ) an additional —i6 prescription has 
been added. As will be shown in the next paragraph, this does not have any effect on the value 
of the four-point function. It does make possible, however, to express the integrals in (^8|) in 
terms of the functions Ii and I2 as in (^l|). 

It is easy to sec, that the limiting procedure as explained above does lead to an unambigu- 
ous result. One might in principle worry that limits — > m| taken from above and below, 
corresponding to the limits a 00 and a — 00 respectively, would lead to different results. 
The question is most conveniently settled if the Dq functions in Eq. ( |50| ) are replaced by the ex- 
pressions given in Eq. (^). Once the limit nii —> m2 is taken, the first factors of the integrands 
have the same limiting value. For a large, thus the leading term is 

/•oc /•! 

Dq ^ —a I 2dX / dy[—2aprP2y + 2av-piX + 2aprp2 — iS]^^ X 
Jo Jo 

[{pi ~ P2fy^ + X^ + 2v {pi - p2)yX + 

+ P2X + {-{pi - P2f +ml- ml)y + - i5]-^ 
1 (53) 

/"OO /"l ^ ' 

+a I 2dX I dy[—2api-p2y + 2av-piX + 2aprp2 — iS\^'^y~ 
Jo Jo 

[ply'' + X'-2v p2yX + 

+ P2A + {-pI + ml - ml)y + - i(5]-\ 



with P2 = 2{v-p2 + A). After collecting the two integrands in (|53| ) the first factor in the 
integrands cancels and one finds 



1 

2„,2 I \2 



Do-^~ 2dX ydy[(pi-p2Yy' + X' + 2v{pi-p2)yX + 
Jo Jo 

+ P2X+{-{pi-p2)'+ml-ml)y + ml-i5]-^ x (54) 
[ply' + X'-2v p2yX + 

+ P2X + {-pI + Jn\ - ml)y + mj - iS]^\ 

This result exhibits clearly the fact that (i) the limit ml ml is independent of whether it is 
taken from above or below and (ii) the limit is independent of the size (or even the sign) of the 
infinitesimal parameter in the first terms of the integrands in ( |53| ) . 

When the momenta pi, p2 satisfy (mi — m2,3)' < Pi 2 < (™i + "^2.3)^, rendering a complex 



a, the procedure outlined above in (pq)-(pd) cannot be applied directly. Starting from (35), we 
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then use the propagator identity (||) on the last two propagators in (|35|), where we set a such 
that P = (pi + a{p2 — Pi))'^ — 0. This has a real solution for a since pi and p2 are timelike as 
has been assumed at the beginning of this paragraph. Changing the notation slightly we then 
have for the scalar four point function (omitting the —iS prescription in the notation) 

^' 

d^q- 



{2n)^J ''{vq~A)[q^~ml][{q+p,)^-ml][{q + l)^-AP] 

1 — a 



d'^q- 



(27r)4 J ^{vq-A) [g2 _ mf] [{q + p^Y - ml] [{q + If - A'P] ' 
with a' the (real) solution of 

{pi+a'{p2~Pi)f = (56) 

and 

I ^ Pi + a'ip2 - pi), (57a) 
= (1 - a')ml + a'ml - a'{l - a'){p2 ~ Pif ■ (57b) 

The integrals in ( |55| ) can now be solved using the procedure outlined above (|3^)-(|5|), once we 
permute the last two propagators with I taking the role of pi in (|35|). Note also, that a' solves 
quadratic equation (do) that in general has two solutions. The final result for the four point 



function Sq does not depend on which of the two solutions is taken in (55). This fact can be 
exploited in the numerical implementation as a useful check. 

The four point function has already been calculated before for a special case of mi = m2 = 
m^, Pi=P2 = and Pi — P2 = Mv^ (see Eq. (All) of [0). It has been checked numerically 
that the two solutions, the one given in this paper and the solution of |l^, agree for this special 
case. A number of other numerical tests have been performed. The direct numerical integration 
of (^) and the evaluation of analytical result given in this paper have been found to agree 
numerically. It has been also checked that the results do not depend on which of the two 
solutions for a, (3 ov a' is taken. The solution for the four-point function calculated in this 
paper has also been checked numerically to have the branch cuts as required by analyticity and 
unitarity. 



6 Conclusions 

In this article the general expressions for the diniensionally regularized scalar two-point, three- 
point and four-point functions with one heavy quark propagator were derived. The calculation 
proceeds along similar lines as in the case of the one-loop scalar functions of the conventional 
perturbation theory. There are, however, also a number of differences that result from the 
appearance of the heavy quark propagator. Since this is linear in the integration variable, 
one is forced to use a modified Feynnian parameterization (||). In the calculation of the four- 
point function this prevents one to use the elegant formalism of jl^,|^6|. A repeated use of 
the propagator identity (|^) then leads to the final expression for the four-point function. The 
simplification in the calculation is, that the heavy quark velocity is a timclikc vector, which 
than leads to real solutions of the quadratic equations (^8|) , ( ^ ) . The complication on the other 
hand is the integration over infinite range that stems from the Feynman parameterization (|^). 

The five-point as well as the higher-point functions can be expressed in terms of the scalar 
functions given in this paper using the standard procedure p5[,]T9|. Consider for instance 
the case of five-point scalar function. This is a function of four vectors, v and pi,P2,P3,- 
The five-point function is first multiplied by v^Cap-yS and then antisymmetrized in all five in- 
dices. The resulting tensor is zero, because there is no antisymmetrical tensor with five in- 
dices in four dimensions. Then the tensor is multiplied first with piaP2/3P3-y<lS and finally with 
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v^^a' p'i'S'Pi P2 P3 ■ Using the decomposition of the product of two Levi-Civita tensors in 
terms of the Kronecker delta functions and expressing the scalar products q-pi in terms of the 
propagators {{q +Pi)'^ — mf^i) and {q^ — ml), the five-point function can be expressed in terms 
of the four point functions. The tensor functions can also be expressed in terms of the scalar 
functions using the algebraic reduction 

In the numerical implementation of the expressions given in this article, further care has 
to be taken regarding the numerical instabilities. Such numerical instabilities can for instance 
arise, if one of the solutions of the quadratic equation is much smaller than its coefficients. There 
is also a possibility of a cancellation between the dilogarithmic functions, when the values of 
the dilogarithms separately are much larger then their sum. These difficulties can be dealt with 
along the lines of Ref. . 
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A Dilogarithm and hypergeometric function 

In this appendix we list some properties of the dilogarithm and the hypergeometric function 
used in the rest of the paper (for other properties consult e.g. ||l8| ). 
The dilogarithm or Spence function is defined as 

Li,(:,) = - /'di Mi^. (58) 
Jo ^ 

The cut for the logarithm along the negative real axis translates into the cut for the dilogarithm 
along the positive real axis for x > I. For x on the positive real axis, a; > 1, dilogarithm is 
calculated using the following prescription Li2(x) — > Li2(a; — ie). Note as well that Li2(0) = 0. 
Useful identities valid also for the complex arguments (not equal to zero) are 

Li2(x) = -Li2(i) -i7r2-l[ln(-x)~27riC(x)]', (59a) 
Li2(a;) = - Li2(l - x) + ^n^ - ln{x) ln(l - x), (59b) 

where 

I ; otherwise. 

The hypergeometric function for complex argument \z\ < 1 is defined in terms of the series 

^0 7---(7 + »^)(»^ + l)! 

with 7 not equal to zero or negative integer. Note that the series terminates if a or /3 are equal 
to negative integer or zero. If either of them is zero then 

2Fi(0,/?;7;z) = 2^^i(a,0;7;^) = l. (62) 

For z outside the unit circle the values of the hypergeometric function can be obtained through 
analytic continuation. We make a cut in the z plane along the real axis from z = 1 to z = 00. 
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Then the series ( plD will yield, in the cut plain, a single valued analytic continuation that can 
be obtained using the following identity (other similar transformation formulas can be found in 
e.g. 



2Fi(a, /3; 7; z) - ^S^i^ (-z)-°2Fi(q, a + 1 - 7; " + 1 - V^) 



r(/3)r(7-a) 
r(7)r(a-/3), 



(63) 



p. (-z)-'^2Fi(/3,/3 + l-7;/3+l-a;l/z). 

1 (aji (7 - (i) 

The integral representations of the hypergeometric function include 



2Fi{iy, ly - fj.;v - ^ + 1; -1/ (/3z)), 



(64a) 
(64b) 



/o (l+/3a;)- ^ 

where Eq. (64a) is valid for Re > Re while Eq. ( |64h| ) is valid for the case, when arg(l+/3z) < 
TT and Re/z > 0. 



In this appendix we will express the integrals appearing in (|j 
First we review the derivations given in jTsl. Consider 



B Reduction to dilogarithms 

P^ ) in terms of the dilogarithms. 

i?(Ai,Ao)= / dA-^[ln(A- Ai)-ln(Ao- Ai)] 
Jo A - Ao 

[InA-ln(Ao-Ai)], 



l-Ai 

d\- ^ 

Ai A — Ao + Ai 



(65) 



where Aq.i may be complex. The residue of the pole of the integrand is zero. The cut of the 
logarithm is along the negative real axis, so for Ai not real, the cut is outside the triangle 0, 
— Ai, 1 — Ai. The integration path can thus be deformed to (for Ai real, this statement is trivial) 



l~Ai 



d\ 



-Ai 



1-Ai 



d\ 



-Ai 



d\. 



Making the substitutions A = (1 — Ai)A' and A = AiA' we obtain 



i?(Ai,Ao) = / d\ 


d\ 



d 



d\ 



In 1 - A 



1-Ai 
Ai — Aq 

Ai 
Ai-Ac 



[lnA(l-Ai)-ln(Ao-Ai) 



)][ln(-AAi)-ln(Ao-Ai)]. 



(66) 



Since A is positive real, none of the arguments of the logarithms crosses the negative real axis. 
After integration per partes 



^(Ai, Ao) 



Lio 



Lio 



Ai-1 
Ai — Ao 

Ai 
Ai — Ao 



In 
In 



1- Aq 
Ai — Ao 

-Aq 
Ai — Ao 



[ln(l-Ai)-ln(Ao-Ai)] 
[In(-Ai) - ln(Ao - Ai)] . 



(67) 



This can be further simplified using ( 59b ) 



i?(Ai,Ao) =Li2(-^V) + ^/f -Ai,-^— ) +2^iRe(-nAo-Ai) 
VAo — Ai/ [ \ Ao — Ai / 

- Li. (^) - kl - A. + 2.. Re(-)(Ao - A,) 



Ao 



In 



Ao — Ai 
Ao-1 



(68) 



Ao — Ai 
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with r] defined in ( pi] ) and Re^~^(a;) defined in (||). Note that this result differs slightly from the 
one in [ p5| as it is defined also for the arguments lying on the negative real axis. The extension 
to negative real arguments was not necessary in Jisf as then the Ao was always real. This is not 
the case in the calculation of the four point function with one heavy quark propagator presented 
in this paper, as the Ai and A2 in (^8|) can have nonzero imaginary parts. The momenta and 
the masses in the calculation can then be chosen such, that one of the arguments appearing in 
( |68| ) can lie on the negative real axis. 
Next we turn to the integral 



S'3(a, 6, c, Ao) 



A — Ao 



[ln(aA^ + 5A + c) - ln(aAo + b\o + c)], 



(69) 



with a real, while b, c, Aq may be complex but such, that the imaginary part of the argument 
of the logarithm does not change sign for x S [0, 1] (also Imc ^ 0). 

Let e and 6 be infinitesimal quantities that have the opposite sign from the imaginary part 
of first and second argument of the logarithm respectively. That is, the signs of the arguments 
are as given by —ie and —id. Then 



^3 = / dX-^ 

lo A — Ao 



[ln(A - Ai)(A - A2) - ln(Ao - Ai)(Ao - A2)] 



77I a — le. 



, — iS 



In 



Ao-1 
Ao 



(70) 



with Ai,2 the solutions of oA^ + bX + c — a(A — Ai)(A — A2). Next we split up the logarithms, use 
the fact that the imaginary part of (A — Ai)(A — A2) has the same sign as the imaginary part of 
c/a and use the definitions of i?(Ai, Ao) (pq) to get 



53(0, 6, c, Ao) =i?(Ai,Ao) +i?(A2,Ao) 

+ 77(-Ai, -A2) - 77(Ao - Ai, Ao - A2) - 7j\a- le, 

with e and S defined before Eq. (|70|). 
For future reference we also define 



iS 



In 



Ao-1 
Ao ^ 



(71) 



£'2(6, c, Ao) 



A — Ao 



[ln(6A + c)-ln(feAo + c)], 



(72) 



with b real and c, Ao possibly complex (Imc ^ 0). Defining as above infinitesimal parameters 
e' and S' to have signs opposite to the imaginary parts of the first and the second argument of 
the logarithms respectively, we obtain 



:(6, c, Ao) = Ri-c/b, Ao) - v{b - -^^) 



In 



Ao-1 
Ao ■ 



(73) 



Next we turn to the integrals appearing in the calculations of the three-point and four-point 
functions with one heavy quark propagator. Consider first 



/i(ai, a2, 61, 62, ^3, ^4,ci, C2, C3, C4, Ao) 



1 r fesA-fcs 
"At ^ In-T"; 

A — Ao L biA + C4 



In 



fliA^ + biX + ci 



a2A2 + 62A + C2 



(74) 



with ai^2 and fei,...,4 real, while Ao, C1...4 may be complex but such that Imcilmc2 > and 
Imc3 Imc4 > 0. Also, the coefficients are such, that for A = Ao the two logarithms are equal, so 
that the residuum of the integrand is equal to zero. Such an integral appears in the calculation of 
the four-point scalar function (^). To reduce the integral Ii to the integrals S2, S3 we add and 
subtract the values of the logarithms at the pole. Since the numerators and the denominators 
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of the logarithms in ( |74| ) have imaginary parts of the same sign, we can spht the logarithms. 
Additional 77 terms appear, however, when we split the logarithms with A set to Xq. As the 
result we get 



h = S'2(&3, C3, Ao) - S2{b4, C4, Ao) - 5*3(01, 61, ci, Ao) + 53(02, 62, C2, Aq) 

^ ^ 2TriRe^-\a2Xl + b2XQ 

Ao-1 



ri(^aiXl + &1A0 + ci, 
- v[h>^o + C3 



64 Ao 



C4 



62 Ao + C2 

- 27riRe^"^(64Ao + C4) 



C2) 



(75) 



In- 



An 



with 77 defined in (0) and Rc^") in (|). 
Next consider the integral 

^2(ai,a2,5i,/i,g2,/2,Ao) 

DO 

rfA- 



A- An 



In 



A-ai , A2+giA + /i 



A — 02 



In 



A2 + ff2A + /2 



(76) 



with (7i_2 real, while Aq, 01^2, /i,2 may be complex with the restriction Imailma2 > 0, 
Im/i Im/2 > 0. Then the logarithms can be split without introducing rj terms, independent of 
the value of A as long as this is real. Also the arguments of the logarithms in ( [76| ) are taken to 
be the same for A = Aq, so that the residuum of the integrand is zero. Such integrals appear in 
the calculation of the three-point function (pG) and in the calculation of the four-point function 



|). We rewrite the integral (76) as 



I2 = 



I?'xh^A 



In _ In (^-^i)(^-^2) 

A - 02 (A — ci)(A — C2) 



with 



A'+5iA + /i = (A-&i)(A-62), 
A'+52A + /2 = (A-ci)(A-c2), 



(77) 



(78) 



where Im bi Im 62 < 0, Im ci Im C2 < 0, Im(&i52) Im(ciC2) > as can be seen from the constraints 
on .91,27 /i,2- Then the logarithms can be split up in the sum of the logarithms with arguments 
linear in A. 

To the integral (|77|) we add logarithms with A set to Aq and then split the logarithms 



- In ~ _ In (^0 ~ ^i)(Ao - ^2) 
Ao - 02 (Ao - ci)(Ao - C2) 



(79) 



where Ki are the coefhcients 01,2, &i,2) ci,2 with p{Ki) — 1 for oi, ci,2 and p{Ki) = —1 for 02, 61,2- 
There is also a sum of 77 terms that we do not write out explicitly, but just denote by 77', as it 
will be reabsorbed in the final result. Note also, that in the case of Ao — Ki real and negative 
the logarithm is calculated using the prescription Xq — Ki ^ Xq — + iS, with 5 a positive 
infinitesimal parameter (see also (P)-(p^)). The integral is then 



^2 = V p{Ki) / 

JO 



dX 



1 /"OO 

[ln(A -K,:)-ln(Ao 7?' / T T' 



(80) 



The separate integrals are divergent so they have to be regulated. We use the cutoff M that is 
sent to infinity at the end of the calculation. Note also, that there is no problem with the pole 
in the last term even if Aq is real, as then 77' is zero. 
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The regulated integrals are then 



dX- 



A- Ar 



■[ln(A - Kt) - ln(Ao - Ki)]. 



(81) 



Let us from here on first assume, that Aq — Ki is not negative real. Changing the variable 
A — A/A' and using the calculation of R{\i, Aq) (|65|), (|68[) we get 













In f A - ^) 




Jo A- If 


V mJ 




Mj 



Li? 



Ao 



Ao 



— Lio 



Ao -M 
Ao — Ki 



Ao 



In- 



-77 1- 



M 



M' Ao - 



In- 



An - M 



(82) 



For M big enough the last term is zero. The M dependent dilogarithm can be transformed 
using relation ( 59a ) 



T- . -M . Ao - 1 
L12 T = - L12 — -TT 

Aq — Ki —M 6 



'-^In^ 



M 



Ao — Kj 



(83) 



The argument of the dilogarithm on the right-hand side goes toward zero as M ^ 00, so that 
in that limit the dilogarithm vanishes. Next we split the logarithm in the last term and write 



In" 



( ^ ) = \n M - 21nAf ln(Ao - k,) + ln^(Ao - K^). 

\ Ao - Ki / 



(84) 



The first term gives zero once summed over in (|80|), while the second term cancels against the 
77' term in (|80|). Leaving the case of Ao — Ki negative real to the reader, the final result is 



I2 =Y^piK,)I^U2 



An 



Ao 

2/ 



l(-Ki,- — )+27riRe' ^(Ao - k,;) 

\ An — Ki / 



In- 



An 



+- In (Ao - Ki)- ln(Ao - Ki) In(-Ao) } , 



(85) 



with 



p(k») = 



1 \ Ki 



Ol, Ci,2, 
= 0,2, 5l.2. 



(86) 



Note that this solution applies also for the case encountered in the calculation of the three point 
function (^), when ai = 02 = 0. Then the terms containing ai.2 cancel each other, so they can 
be dropped altogether for the case of Eq. (^0|). 

There is one more point worth mentioning regarding the expression (|85|). One might think 
that problems could arise for Ao — Ki negative real or Ao/(Ao — Ki) real as then one has to 
deal with the cuts in the logarithm and the dilogarithmic function^. We use the prescription 
for the arguments lying exactly on the cuts of the functions as described before Eq. (||) and 
after Eq. (|58|). One could as well use a different prescription, with infinitesimal parameter e in 
(^), (HI) taken to be negative, and with appropriately adjusted rj and Re^~' functions. It has 
been checked numerically that the result ( ^sf ) does not depend on which prescription is used. 
Thus the result (|8|) is valid for any complex Ao , Ki independent of the prescription used for the 
arguments lying on the cut. 



'Note that there exists such a combination of parameters v, pi.2, A and mi, 2, 3 in (M) that Aq — Ki in 



negative real for some i, as can be seen from definition of ai, . . . ,d2 (^), definition of Ai, 
for the four-point function ( ^ ) 



and the expression 
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For the special case of Ao real the result ( p5D simplifies considerably. The rj term is then 
zero. Also the last term in (|85|), that arises from the ij' term in (|80|), then sums up to zero. For 
Ao real we have 



p{Ki) 



Lio 



Ao 



An 



- ln^(Ao - Hi) 



(87) 



with Ki and p{Ki) as in (p6[). 

Of special interest is the case of Ao almost real, i.e. 



An 



An^ + 1(5', where An'^ is the real 



part of Ao and S' an infinitesimal (not necessarily positive) parameter. One can of course still 
use the solution (|85|). The problem is, however, that for both Ki and Ao almost real one has to 
keep track of the relative sizes of the infinitesimal imaginary parts. This complication can be 
avoided by the following procedure. First we set Aq in the second line of ( |79| ) equal to its real 
part. By doing this, the arguments of the logarithms can cross the negative real axis, which is 
compensated by a new sum of 77 functions, 77". Then instead of (pfl) we have 



E 



p(kj) / d\ 



1 



A- A?; 



iS' 



;[ln(A-«,)-MA5'=-«.)]+?7" 



dX 



A-Ao 



In the first integral d' can be safely put equal to zero as the resulting integrand has vanishing 
residuum, with the logarithm in the numerator being an analytic function in some neighbourhood 
of the pole (since Im(Ki) ^ 0). The integral thus docs not depend on how we avoid the pole (i.e. 
6' can be positive, negative or zero). In the second integral one has to keep the imaginary part 
of Ao. 



The final result for almost real Aq = Afj^ + i5' is then 



Li2 



^ln'(A5^ 



K,) - ln(Ao'' - Ki) In(-Ao) 



(89) 



where in the last logarithm Ao is kept together with its infinitesimal imaginary part. 
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